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TURBULENT BOUNDARY LAYER ON A YAWED CONE IN A SUPERSONIC STREAM 

By Willis H. Braun 


SUMMARY 

The momentum integral equations are derived for the boundary layer 
on an arbitrary curved surface, using a streamline coordinate system. 
Computations of the turbulent boundary layer on a slightly yawed cone are 
made for a Prandtl number 0.70, wall to free-stream temperature ratios 
of l/2, 1, and 2, and Mach numbers from 1 to 4. Deflection of the fluid 
in the boundary layer from outer stream direction, local friction coeffi- 
cient, displacement surface, lift coefficient, and pitching- moment coef- 
ficient are presented. 


INTRODUCTION 

The flow of a laminar boundary layer over a cone in a supersonic 
stream has been well established, not only when the cone is alined with 
the free stream (refs. 1 and 2) but also for certain perturbed motions. 
Moore (refs. 3 and 4) has computed the boundary layer on a cone at angle 
of attack, and Illingworth (ref. 5) has found the flow when the cone is 
spinning. Recently, the laminar boundary layer for a cone that is both 
spinning and at angle of attack has been investigated (refs. 6 and 7). 

The corresponding treatment of the turbulent boundary layer on a 
cone has only been begun. Van Driest (ref. 8) and Gazley (ref. 9) have 
studied the cone at zero angle of attack. They have found transforma- 
tions relating boundary layers on cones to those on flat plates. 

Van Driest* s method represents the turbulent stresses according to the 
mixing-length theory, while Gazley uses the momentum- integral equations. 

The investigation of the cone at yaw (attack) or in spin is impeded 
by the lack of any generalization of the turbulent stress representations 
to three-dimensional boundary layers. Moreover, it appears from experi- 
mental measurement (see ref. 10 for review of experimental results) that 
the "principle of independence** may not hold for three-dimensional tur- 
bulent boundary layers. That is, on a yawed cylindrical surface, the 
flow in the plane normal to the axis may not develop independently of 
that along the generators, as in laminar, incompressible flow. From 
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these considerations it appears that a resolution of the equations of 
motion along a body coordinate system will not, in general, lead to a 
tractable problem. More appropriate for this case is the integral method 
as applied by Mager (ref. ll) , which follows the development of the 
boundary layer along a streamline, assuming that the thickness, skin 
friction, and other properties behave nearly as in plane flow. 

The following development extends Mager ! s use of the integral method 
to compressible flow over arbitrary curved surfaces. The specific case 
of the turbulent boundary layer on a cone at small angle of yaw is con- 
sidered. The forces exerted by the boundary layer on the cone and the 
heat transfer to the surface are computed. 


SYMBOLS 


B 


A 


defined in equation (19) 
defined in equation (18) 


C 


constant 



frictional drag coefficient 

local friction coefficient (defined in eq. (l2)) 


pitching-moment coefficient 
defined in equation (51 ) 
defined in equation (42) 


frictional lift coefficient 


D 


F 



frictional drag force 



frictional lift force 


g 


determinant of metric coefficients 



metric coefficients 


I 


integral 


K 


constant 


L 0 ,Li,Ni,Pi,Ro,Ri 


defined in equation (40) 
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M 

Mach number on surface of unyawed cone 

M s 

Mach number on surface of yawed cone 

Mi 

Mach number ahead of shock 


pitching moment 

n 

index , equation (15) 

Pr 

Prandtl number 

P 

pressure 

r 

distance from apex of cone 

St 

Stanton number 

s(i) 

arc length in i*'* 1 direction 

T 

temperature 

% 

t 

time 

* u 

velocity component in direction of cone ray 

uVco 

velocity component in i th direction 

w 

circumferential velocity component 

X, z 

constants defining potential flow 

X 1 

coordinate in streamline system of coordinates 

a 

deflection angle at wall 

a^ag 

defined in equation (45) 

r 

ratio of specific heats 

A 

displacement surface height 

& 

boundary- layer thickness 

5 1> B 2 

displacement thicknesses, equation (ll) 

m 

e 

yaw angle 
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n 

0 

d ll’ 9 21> 9 22 

X 

V 

K 

p 

<P 

* 

Subscripts : 
m 
s 
w 

Superscript: 


parameter defining pressure in potential flow 
cone half- angle 

momentum thicknesses, equation (ll) 

yaw parameter 

viscosity 

kinematic viscosity 

parameter defining density in potential flow 

density 

viscous stress 

angle about wind axis 

angle about body axis 

mean value 
stream value 
wall value 

value on unyawed cone 


INTEGRAL EQUATIONS ON A CURVED SURFACE 



x° 

2 u 2 / 

L 

/ J 

// 


u c 


u 


a 


/ 


y 


(a) 
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Consider the boundary layer of thickness 6 on the surface S of 
sketch (a). The projection of the stream velocity u s on the surface 

determines the x 1 direction, (it is assumed that u s is the velocity 
a potential flow would have at the surface.) The direction in the plane 
of the surface normal to x 1 is designated x 2 . With the choice of x° 
normal to the surface the result is a right-handed orthogonal system. 
There can, in general, be a pressure gradient in the x 2 direction. The 
effect of the pressure gradient is much stronger upon the low-energy 
fluid in the boundary layer than upon the outer stream. Consequently, 

u does not vanish in the boundary layer, and the turning of the stream- 
lines within the boundary layer is greater than at x° = 5. The angle 
of deflection between the outer streamline and the limiting direction of 
the streamlines at the surface is denoted by a. 


and g 22 
By choosing 

§00 = 1 * 


The element of arc ds in the x- system is given by 
(ds) 2 = go 0 (dx°) 2 + gqqCdx 1 ) 2 + g 22 (dx 2 ) 2 
where the g^ are the metric coefficients. Of these, g^ 


are functions of 


*°, 


x , x^ determined by the surface S. 


x° to be the normal distance from S, it is correct to set 
With the definition 


= g 


3 iJl 


= g ll g 22 


the equation of continuity becomes (ref. 12) 



-A_ (pu 1 ) + i pu 1 

Ux* 2 





0 


(1) 


The momentum equations are (as derived from ref. 12) 



1 dp 
®ii dx^- 



i = 0,1,2 


( 2 ) 


In equation (2) u 1 , is the covariant derivative of the contravariant 

vector u 1 with respect to x L , and is the stress tensor excluding 

pressure. 
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It is more convenient later to have (l) and (2) in terms of the 
physical components of the velocities and stresses. The physical com- 
ponents are given by 

u (i) = V^ii u± (3) 



(4) 


and the arc lengths in the three coordinate directions are 


ds(i) » V§ii to 1 


(5) 


If, in addition, the results of the appendix are used to choose the 
dominant turbulent force terms, equations (l) and (2) become, respectively, 


t 

i=0 L 


^ + -nr ln V5 




§s( 


= 0 


and 


£ Itfry] + pu(0)u(1) * + 

i=0 L * J 

- pu2(2) Sifiy ln ^22 = 


[ pu (b + pu (°) u ( 2 ) - . s(0 ) 22 + P u d) u (2) — 

2 ,^ * ln ViIl d<02) dp 
■ pu (1) ■ -ssrzr" = ^sXoT-" 

Outside the boundary layer, equations (7) and (8) become 


Ps u s 


du g 

^iy 


dp 

'ssm 


(6) 


(?) 


( 8 ) 


(9) 


and 

2 ^ ln VSn d P 

p s u s - ds(2) = ■d^2T 


( 10 ) 
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Equations (7) and (8) are now to be integrated across the boundary layer 
making use of (6), (9), (10), and the following definitions of displace- 
ment and momentum thicknesses: 



The result is 


^11 


+ e 


11 




ln Ps u s 



a 

c>s(l) 


ln u s 


ae 


21 


^2j 



as* 

^2) 


- 6 


* £ 

2 as (2) 


ln P s u s 



T ( Q1 )w 

Ps U l 



( 11 ) 


( 12 ) 


ae 


” 21 , n O Z , 66 t ~ On 2 / 

+ 0 21 ln Ps u s§22 + + 0 22 ^ P s u g ->/gll 

a 


ae 


22 


(e 


11 


+ 5*) 


a s ( 2 ) 


1 T (02) w -L 

ln V®11 “ 2 Cf t(01) w " 2 Cf tan a ( 13 ) 


In deriving (12) and (13) use has been made of the fact that the metric 
coefficients change very little over the thickness 5 of the boundary 
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1 +CT ° 


layer, so that (ref. 12) 

sifn ln ^ ' s^iy ln V5l| x0=0 [ 3 

SsfaJ in VS^ = 5^7 in [l +e (|)' 



±ln Vi ~ - 
Sx° R 


(14) 


where in each of equations (14) R is some radius of curvature. Terms of 
order 5 /r compared with unity have been neglected in (13). 


Equations (12) and (13) govern the development of a compressible 
boundary layer on a curved surface on those portions of the surface where 
the radius of curvature is very much greater than the boundary- layer 
thickness. In this respect the foregoing derivation resembles that of 
Timman (ref. 13) for laminar, incompressible boundary layers in con- 
trast to that of Mager (ref. 14) for incompressible boundary layers over 
axisymmetric bodies at yaw. 


VELOCITY PROFILES AND FRICTION COEFFICIENT 

The five boundary- layer thicknesses (ll) that arise in the equations 
of motion (12) and (13) are functions of the two velocity profiles in 
the boundary layer u(l)/u g and u(2)/u s , and of the density ratio 
p/p s . These three unknown variables along with the friction coefficient 
Cf and the deflection angle a give a total of five. Therefore, hy- 
potheses about the boundary- layer characteristics must be added to the 
two equations of motion for the system to be solvable. The only hypoth- 
eses available are those suggested by studies of the boundary layers on 
flat surfaces where the streamwise pressure gradient is negligible. It 
is convenient to write these hypotheses here and to discuss in a sub- 
sequent section whether they are appropriate to a cone at yaw. 

With respect to the velocity profile in the direction of the free 
stream, the customary assumption of a power law is madej as in planar 
flow, 

^ - (€ /n 

The index n is known to vary roughly as the logarithm of the Reynolds 
number (ref. 15), although to simplify computations it is usually kept 
at a constant value of 7 to 9. 
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The component of crossflow in the boundary layer u(2)/u s must 
vanish both at the wall and at the edge of the boundary layer , as shown 
in sketch (a). An expression which yields a qualitatively correct pro- 
file and which is also quantitatively accurate for slow flows (ref. ll) 
is 


u(2) 

u s 



tan a 


(16) 


For the density ratio p/p s there is available under the condition 
of negligible pressure gradient the quadratic in the velocity profile 
that holds in plane flows (e.g., ref. 8): 


where 



(17) 

(18) 
(19) 


The subscript s refers to conditions outside the boundary layer, and 
T y is the (constant) temperature of the wall. The factor Pr p /^ is an 
empirical correction (ref. 16). 

Again, using a property of planar flow, the friction coefficient is 
replaced by a generalized Blasius formula 



The subscript m refers to reference conditions in the boundary layer. 
Eckert (ref. 17) gives for a reference temperature 

T m - i (Tw + T s ) + 0.22 Pr 1 / 3 3L -£- i T s (2l) 

The value of K is generally taken as about 0.045 for a flat plate and 
— — (0.045) for a cone alined with the stream. 

V3 

The relations (15) through (2l) reduce the number of dependent vari- 
ables to two, the thickness 6 and the deflection angle a. Stream 
Mach number and ratio of wall to stream tenperature are the parameters. 
Further progress in the solution of the equations of motion (12) and (13) 
requires knowledge of the stream of a given conf iguration (i.e., Mg). 
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STREAM COORDINATE SYSTEM ON A SLIGHTLY YAWED CONE 


The coordinate system on the yawed cone is composed of the 
potential-flow streamlines on the surface and their orthogonal trajec- 
tories, A comparison is made in figure 1 of the flows on the surface of 
an unyawed cone and a yawed cone. Whereas the surface streamlines coin- 
cide with the elements of an unyawed cone, on a yawed cone they are de- 
flected by the circumferential pressure gradient. 


The potential flow may be found for small yaw angles e by use of 
the tables of reference 18. (The definition of e used here and in 
ref. 18 is the negative of that of some treatments.) 


In terms of the angle cp about the wind axis, and to first order 
in e, the radial and circumferential components of velocity are given 

by 

u =s u + ex cos cp (22) 

w = ez sin cp (23) 


where x,z,u are coefficients tabulated in reference 18 against cone 
angle and Mach number. The velocity u is recognized as the radial 
velocity at no yaw (s = 0). When making computations it is necessary 
(ref. 3) to replace the values of z appearing in reference 18 by 
z = -z (tables) -2x/sin 0. At the surface of the cone it is allowable 
to first order in e to replace the angle cp by the cone angle \|r 
(ref. 19). Then the path of the streamlines on the surface of the cone 
is given by 

dr 

u _ dt 
w ~ r sin 0 

dt 



Integration of (24) yields 


f 



1 - cos 

1 + COS l|f 



- 2L sin 


0 


" 2 sin 0 


r u (sin \|r) u 


(25) 


2 

where each value of the parameter x is associated with a streamline. 
The exponent e/sin 0 appearing in (25) is the parameter describing the 
perturbation of the system from zero yaw, and is hereinafter denoted as 

X: 


X = 


e 


sin 0 


(26) 


4372 


CP-2 back 4372 


NACA TN 4208 


11 


Let the function f (x 2 ) be chosen so that at zero yaw (X = o) x 2 = -cos \|r. 
Then 

f(* 2 ) - ^4 (27) 

1 - X* 

Equations (25) and (27) taken together are an implicit transforma- 
tion from r ^ \|r to x . Their combination yields 

x 2 = -cos * - sin 2 * In r z A* (sin 0 (28) 

Equation (28) describes the streamlines on a yawed cone to first order 
in the yaw parameter X. The new coordinate x 2 runs from -1 at the 
top of the cone to +1 at the bottom. 

The curves orthogonal to the streamlines are defined by 

dr ez sin * 

r sin 0 d* u + ex cos * 

which upon integration yields 

x^ = r(l - X sin 2 0-^-cos *) (29) 

On the left of (29), x^ is a parameter associated with the family of 
normal curves or, alternatively, a coordinate measured along the stream- 
lines. Equations (28) and (29) are the transformation from the unyawed 
to the yawed coordinate system. 

The inverse transformation is readily found to be 

r a x-*-(l - X sin 2 ®-=-x 2 ) (30) 

cos * = -x 2 - X[l - (x 2 ) 2 ]{^ In x 1 - -- ^_ n 9 In |l - (x 2 ) 2 j| (3l) 

On the surface, 

(ds) 2 a (dr) 2 + (r sin ® d*)^ 

= SuOax 1 ) 2 + §22 O^x 2 )^ (32) 
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Using (30) and (3l) in (32), it is found that 

= 1 - 2X sin^0-§- x 2 


u 


*22 


= (* lsln e > 2 (l - 2X*2(4- 

1 - (x 2 ) 2 \ l u 


sin 2 © + 4 In x 1 - — sin 0 


u 


X sin 0 
u 2 


In [l - (x 2 ) 2 ]}) 


(33) 


(34) 


The density and velocity at the edge of the boundary layer are also 
required in the equation of motion. They are obtained from reference 18 
to first order in X as follows: 


P s = p(l 


T s = T 


1 + X sin 


and 


M 2 = M 2 


1 + X sin 01 


. - 2X X sin 0 x 2 ) 
u 

(35) 

K sin © „2) 
P 

(36) 

T) sin 0 x 2\ 
P 

(37) 

that 


<WH 

(38) 

VP u pJ J 

(39) 


In the preceding equations, a bar over a quantity indicates its value at 
zero yaw. 

With the aid of expressions developed, it is possible to show that, 
for a cone at small angle of yaw, the forms hypothesized for the velocity 
and density profiles in the preceding section are indeed appropriate. 

The substitution of (16), (33), (34), and (37) into the equation of 
motion (7) gives, upon neglecting terms of second order in a and X 
and the term of order 5, 


to h (i) 


dx(oi) 

= 
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The pressure gradient is seen to be negligible, so that the velocity pro- 
files (15) and (16) are good approximations. 

Moreover, replacing u(l) by temperature and t( 01) by heat flux 
oyer heat capacity in the preceding equation results in the energy equa- 
tion. This is the condition necessary for Reynolds' analogy between 
stress and heat flux to hold. Consequently, the density (or temperature) 
profile (17), which is based on Reynolds' analogy, can be properly used 
for the cone at small angle of yaw. 


SOLUTION OF THE MOMENTUM INTEGRAL EQUATIONS 

The momentum integral equations (12) and (13) are now to be written 
for a yawed cone at constant temperature in terms of the boundary- layer 
thickness 6 and the deflection angle a. In order to do this it is 
necessary to evaluate first the integrals (ll) using (15) through (21). 
It is convenient to specialize at once to n = 7 in the profiles (15) 
and (16) and to make the following definitions : 


m 


* r 1 m 

Jo i + Bt, - a 2 (; 2 


tm ’ 2 R/f 


Cdt 

(1 + B£ - A 2 £ 2 ) 2 


4A 2 + B 2 


2A ^2 + 2A2 ^ m " 1 ) I m " BmI m-l 


1 + B - A 


V 2 S PrV3 m 2 X 


T. 


w 


B (1 + - 1 


Bi-f [1 Pr 1 / 3 (r- i) m 2 1 

^•wLp P U J 
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Then to first order in X and a, the displacement and momentum thick- 
nesses (ll) become 


c ■#- _ 

°1 T T 

-£-=l-7-=r-I 7 +X sin 0 =- 




= L~ + X sin 0 L-, x 6 


L w 


/ T! £\ x -2 

(p ‘ p) I 7 + ¥8,2 + 2 u A I 9.: 


4 - = a 7 (l 7 - 2I 14 + I 2 i) 


T, 


w 


s a N]_ 


-|i= a^- 7(I 8 - 2I 15 + I 22 ) 


>11 


= a Pi 


7 Y" (I 7 - I 8 ) + x sin 0 7 | 


w 


w 


ft ■ f) 


(I 7 - Is) 


- 2 , 


" B l( I 8,2 ~ T 9,2) " 2 f" A ^S; 2 " I 10,2) 

= Rq + X sin 0 R^_x 2 


(40) 


0 22 2 ^ 

— = flr(a^) 

The integrals I ? , I Q , and (l y - Ig) appearing in (40) are listed 
in the tables of references 20 and 21. The integrals Ig and I10, 
which are required to compute Ig^ 2 B 10,2’ ma y Be obtained from 

the recurrence relation 

B m+1 = =2^ BI m + Im_1 ” m^ 

The combinations of integrals (I7 - 21^ + I 2 p) an< ^ (^8 ” 2B 15 + B 22^ 
should be computed separately, rather than from the recurrence relation, 
to avoid cumulative error. 


4372 


MCA TW 4208 


15 


The Blasius formula (20) may likewise he expressed in terms of stream 
parameters at zero yaw. If a linear viscosity law is used, 

-H- = C — 

>*s T s 


then (from (20)) 



Introducing the thicknesses (40) and the friction coefficient (4l) 
into the momentum integral equations (12) and (13) gives, to first order 
in a and X, 


~T 5 ( R o + ^ sin 0 R i x2 ) + 8 ( R 0 + sin 0 R l x2 ) “T (l - X4- x 2 ) 
ox‘ L x x u 


" (x2 j 2 [l + Ni6a - N^Ba X)| 

sin 0 x 1 Lox^ 

^ ~ ^ x2 j 2 [i + p i 5a - Pi&cc <*(x)l 

sin 0 x x J l_dx^ 


“ K (f) / fe) / t 1 + x 8ln 8(p - t sln s)x2 ] 


and 

8 


— P l5 a + PiSa ± [l + (xj) 

- s£Rq + Lq + X sin 0(R^ + L-^x^J X sin 0 — ~ j x ‘ ^ 

’ aK ft) 1/2 &) 1/4[l +MX)] 


(43) 


(44) 
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In order to solve (43) and (44), the dependent variables a 
are expanded in powers of (small) yaw parameter X: 

o 

tan a = a^X + ctgX + . . . 

8 = 8q + 8-j^X + . . . 

Equation (43) yields equations of zero and first order in X, 
(44) yields a first-order equation as follows: 



, Jl - (x 2 ) 2 „ 1 " N 1 

+ xi sin 0 80 [ ~ ~ 



cbc 1 


^0°! + ^ 


5 0 a i 



sin 0 


(Lq + Rq) 5 0 


Jl - (x ^) 2 

1 



Equation (46) is simply the integral equation for the unyawed cone 
solution of (46) satisfying the condition &(x-*- = o) = 0 is 


u6 0 


v 



Substitution of (49) into (48) yields 


a l = 


sin 0 


z_ 

u 


L 0 + R 0 \ r 

14P ± + 9Eq ) V 


(x R ) R 


and 5 

(45) 
and 

(46) 


(47) 


(48) 
. The 

(49) 


(50) 
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Then using (49) and (50) , the solution of (47) is found to be 




(51) 


Equations (50) and (51 ) show that the boundary- layer thickness has 
its greatest changes at the windward and leeward rays of the cone and 
that the fluid in the boundary layer is deflected most in the horizontal 
plane (fig. l) . 


BOUNDARY - LAYER EFFECTS 

From the results of the previous section it is possible to find the 
effect of the boundary layer on the outer stream and the forces it exerts 
on the body. 


Displacement Surface 


Let A(x^,x^) designate the distance, in the normal direction from 
the body surface, which the potential flow is displaced by the boundary 
layer. Moore (ref. 22) has shown that A is determined by 


div 



(p s u s - pu) ds(0) 


0 


(52) 


where the arrows designate vector quantities. In terms of the previously 
defined boundary- layer thicknesses, (52) becomes 


P S U S (A - S*) = ^2 Ps u s V5ii &2 


(53) 


Equation (53) and the boundary condition at the tip show that, on the 
unyawed cone, the surface A has the value of the usual displacement 
thickness 8^. Substitution of (33) through (36), (40) , (49), (50), and 
(51) into (53) yields 


A 

L 0 5 0 


1 



sin 0 


25 z_ % L 0 + R 0 \ 

9 u *^*1 + 9*0/ 


(54) 
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For a positive value of the yaw parameter *•, A is nearly a circular 
cone whose axis is inclined downward from that of the body. 


Local Friction Coefficient 

Because Reynolds' analogy shows that the local friction coefficient 
is proportional to the local heat-transfer coefficient, it is desirable 
to know the local distribution of stresses on the surface of the cone as 
well as the net forces. Combining (29), (4l), (45), (49), and (5l) 
yields the friction coefficient: 



Using as a measure of heat transfer the Stanton number, 


St = — t m »— 

c p( T wo “ T w)P u 

where q w is the heat transferred at the wall, Cp is the heat capacity, 

and T is the insulated wall temperature, then 
WO 



Frictional Drag Force 

For the net forces on the cone it is convenient to transform the 
stresses from the streamline system of coordinates (x- system) to the 
body system (r,\|r,y, fig. l) • The transformations required are 

T(yr) = t(oi)(|t- Viii + H It - V&22) "1 


T(y*) = r S in 9 T ( 01 )(lr V^I + S. W ^22) 


J 


( 56 ) 


4372 


CP-3 back 4372 


NACA TN 4208 


19 


from which 


x(yr) = x(0l)[l + <?(X 2 )] 


(57) 


R o + L o 


x(yitf) = \t( 01) sin t =- sin ©fl + 5 + g - R 


J 


The force on an element of area da^ on the cone surface is 


= T ij daJ 


The net force in the direction of any unit vector v is 


F(v) = x„v d da^ 


(58) 


If v is in the direction of the cone axis, it has no \|r- component , and 
the drag force is 


F-q = 2 cos 0 sin 0 


IT 


x(yr)r dr d\|r 


The viscous stress appearing in the integral is, by (12) and (57), 
x (yr) = 'F(Ol) /l - X cos \|/ £sin ® F - ^ D - -g sin 2 0 -= 

- - •(* # ♦ f )] } 

The frictional drag coefficient is then 


(59) 


— — p p 

pu^jtr^sm 0 


■ f (cos 0)K (|) 1/2 [^(|) 1/2 ( f)]" 1/5[l +<nx2)] (60) 

To first order, there is no change in the drag coefficient with yaw. 

Comparison with the pressure-drag coefficients of reference 18 re- 
veals that the pressure forces are of the order of a hundred times the 
frictional forces. 
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Frictional Lift Force 

Consider the force in the direction of unit vector v in the plane 
of yaw and perpendicular to the axis. 


The components of v are 



( v r> v i|r> v y) “ ( sin ® cos tj-sin i|/,cos 0 cos i|r) (6l) 

Combination of (57), (58), (60), and (6l) yields for the ratio of the 
frictional lift force to the frictional drag force 




sin 0 - — D 


1 

5 


• ? z 
sin” 1 © — - 
u 


sin 


+ f ( 1 + 5 



(62) 


Frictional Pitching Moment 

Consider the moment of forces about an axis through the cone apex, 
normal to the plane of yaw. The stress T(ry) acts through the apex and 



so cannot contribute to the moment. The stress 'c(i)ry) has a moment arm 
r cos 0, and its component in the upward direction is 

-sin \| i T(\|ry) 


The pitching moment is then 


J( = -2 sin 0 cos 



sin tyT(i|ry)r 2 dr di|r 
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Use of equations (57) and (60) results in 

_ ■£. 

- Fj-,r cos 0 


= -X — tan 0 4(1 + 5 

28 u l 


for the frictional pitching moment. 


R 0 + L 0 \ 
14Pi + 9 Rq J 


(63) 


NUMERICAL EXAMPLES 


Computations of the boundary layer on the yawed cone have been made 
over a free- stream Mach number range of 1 to 4 to show the effect of cone 
angle and temperature ratio. The Prandtl number was taken as 0.7. 


The computational procedure is as follows. For a given Mach number 

and cone angle, the parameters u, x, z, Kfp> and iq/p are obtained 
from reference 18. For a given Mach number and temperature ratio, the 
integrals Iy and Ig are obtained from reference 20 or 21. The in- 
tegrals Ig and I 10 are obtained from Iy and Ig by the recursion 
formula. The combinations of integrals (iy - 2I 14 + Igp) and 
(I 8 - 2I 15 + I 22 ) are functions of M (zero yaw Mach number on the cone 
surface) and the temperature ratio t/t w . They are plotted in figure 2, 

which shows that they are nearly linear in M^. From these parameters 
and integrals, the quantities Lq, L-l, Rq, R-^ P-j_, and N-]_ (eq. (40) ) 

may be obtained, and with them the coefficients of the previous section. 


In figure 3 is shown the parameter (from eq. (54)) 


1 4 ^ i L 1 25 z N 1 L 0 + **0 

14P 1 + 9P 0 


which determines the deflection of the displacement surface axis from 
the cone axis. Figure 3(a) shows the effect of temperature ratio on a 
20° cone, and figure 3(b) shows the effect of cone angle at constant tem- 
perature ratio. A comparison with the displacement surface of a laminar 
adiabatic boundary layer (ref. 3) is made. The parameter is of the same 
order of magnitude for the two boundary layers over the Mach number range 
considered, although slightly higher for the laminar flow. The displace- 
ment surface deflection is slightly greater for low wall temperatures 
because of the increased density in the boundary layer. 
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The deflection of the inner streamline is determined by the param- 
eter (eq. (50)) 


a l 

sin \|r 


sin 



f L 0 + R 0 \ 

Vl 4P l + 9R 0 / 


Figure 4 shows that the deflection is considerably less in the turbulent 
boundary layer than in the laminar. The fluid in the laminar layer, 
having less kinetic energy, is influenced more by the transverse pressure 
gradient. Figure 4(a) shows the effect of temperature ratio, and figure 
4(b) shows the effect of cone angle. When the wall temperature is high, 
the density in the boundary layer is lowered, and greater deflection 
results . 

The change in local skin- friction coefficient ^(0l)/p s u| with yaw 
is determined by (eq. (35)) 


1 

C f cos \|r 


SCfl 

^Jx=o 


F sin 0 - 7 D - i sin 2 © -4 
9 5 u 


Figure 5, a plot of this quantity, shows that the effect of temperature 
ratio is small, whereas the effect of cone angle is appreciable. The 
comparison with the laminar flow again reveals that the turbulent bound- 
ary layer is less affected by yaw. 

The viscous lift coefficient (eq. (62)) is shown in figure 6 . There 
are two types of forces contributing to this coefficient. The viscous 
shear is greater on the top of the yawed cone than on the bottom, re- 
sulting in a lift force in the direction of yaw. The circumferential 
forces always give a lift force opposing yaw. For small-angle cones the 
drag forces grow relatively larger with Mach number until finally the 
resultant lift is in the direction of yaw. 

The pitching moment is shown in figure 7. It is proportional to 
that part of the lift caused by the circumferential forces. It is 
slightly higher for the higher ratios T w /T because of the greater de- 
flection of the low-density fluid in the boundary layer. 


CONCLUDING REMARKS 

The over- all picture given by the foregoing computations is that the 
turbulent boundary layer on the yawed cone behaves qualitatively like the 
laminar layer, but the relative magnitude of the effect is always less. 
Thus, the deflection of the fluid in the boundary layer, the frictional 
lift, and the frictional pitching are relatively smaller because of the 
higher shear stresses. 
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The assumptions entering into the formulation of the problem should 
be kept in mind in interpreting the results. The assumption that the 
skin- friction law for plane flow holds for flows with slightly curved 
streamlines should not be far from the true picture. The choice of a 
power- law profile for the velocity component in direction of the outer 
stream is considered reasonable because the pressure gradients are small. 
But the selection of the proper profile for the transverse component is 
less certain and at the same time critical, for it is this profile that 
determines the three-dimensional effects. The results presented here 
give magnitudes and general trends, but are not necessarily quantitatively 
accurate. 

The analysis gives no information about the separation of the bound- 
ary layer. It is limited to small perturbations, from zero yaw, of the 
potential stream; consequently, there are no strong pressure gradients 
to induce separation. Any form factor constructed from ratios of the 
boundary- layer thicknesses varies at most to order X over the cone sur- 
face, whereas separation is associated with the growth of a form factor 
to some large critical value. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 24, 1957 
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APPENDIX - SELECTION OF DOMINANT STRESS TERMS 

The turbulent stresses arise in equation (2) in the covariant de- 
rivative In a plane flow the normal stresses T fy and the 

shear stresses t ^ t q are too™ from experimental measurement to he of 

the same order of magnitude. In three-dimensional boundary layers, ad- 

p o n l 

ditional stresses T o* T l' T 2' T 2 ar i se that may he smaller, hut certainly 
no larger, than the normal stresses. 

In the x^- equation of motion the stress terms are 


-i- x 1 + JL t 1 + -A- t 1 

§00 0,0 g ll 1,1 §22 2,2 


E 4© + to - to) ♦ + to - to. 


i .(H 


S 22 V dx 2 


+ *12*4 " I 22 T i 


(Al) 


where it is understood that i is summed from 0 to 2. The f are 

mn 

the Christoffel symbols: 


T l = 0 
mn 


for 1 / m / n £ l 


Also, 


r 1 = r 2 =0 

00 00 


and 


§~~ = 1 
00 


Using (A2) and (A3), the right side of (Al) becomes 

+ j + - 

11 dx 1 g 22 Sx £ 


1 1 dT 2 * -fi 1 I 1 + T L fr 1 - r° - JL 

i? ^ 01 °V 10 «li ru ®22 


— r9 2 
g 22 22 


(A2) 


(A3) 


-^ 4 r 22 *i(-i^ r U*i^ 


1 r 1 


JL r 2 

12 g 22 22 


2 

+T 1 


1 

§11 



+ TTc 


1 I 1 

1 O 


§22 


22 


(A4) 
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Consider first the partial derivatives in (A4) . They may be written 
in terms of the physical components and arc lengths, 


(A5) 


1 

’ ^r(io) Sr(llj 

1 8xCl2)l 

V§n 

ds(o) ds(l, 

1 ds(2)J 


In a boundary layer, the stresses vary rapidly in the direction normal 
to the surface, but slowly in directions parallel to it. Thus, 

d 1 d , 3 1 


^iToT 


~ — , while 
o 




and 


ds(2) 


~As ; where As is a distance, 


measured on the surface, over which a stress varies an appreciable part 
of its total magnitude. Then, 


dT(lQ) / dx(n) dx(l0) / cfc(l2) As 
ds(0 ) j ds(l ) ’ ds(0 ) j ds(2) 5 


» 1 


(A6) 


The Christoffel symbols are composed of terms of the form 

1 ^®mm c) 

Therefore, to take a typical term as an example, 

d In g lx 




r 1 ~ 


1 ®n 21 Vi! 


t(21) 


T (2l) 


bs(z) 


VsT As 


(A7) 


=11 


Consequently, 


v& 


11 


8<io 
ds(0 
.2 


g- 


11 


r 1 ~ 

21 5 1 


Thus all terms not involving derivatives of the stresses may he neglected. 
(The fact that some of the derivatives in the Christoffel symbols are 
taken with respect to x° does not make them large, for the g^ are 
nearly constant through the boundary layer.) 

It is seen that the only stress term which remains in the x^ - equa- 
tion of motion is 


1 3r(l 0) 

Vi! ^toT 
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A similar argument shows that in the x equation of motion the only 
term of importance is 


1 dT(20) 

~^(o) 
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Figure 1. - Coordinate systems on unyawed and yawed cones. 
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1 2 3 4 5 

Mach number ahead of shock, Mj_ 


(a) Cone half -angle, 20°. 

Figure 3. - Displacement surface parameter. 
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Mach number ahead of shock, M^ 


(b) Temperature ratio, T w /T, 1. 

Figure 3. - Concluded. Displacement surface 
parameter . 


sin 


32 


NACA TN 4208 



Mach number ahead of shock, M^_ 


(a) Cone half -angle, 20°. 

Figure 4. - Deflection of limiting streamline 
at surface with respect to potential- flow 
streamlines . 
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Mach number ahead of shock, 

(b) Temperature ratio, !• 

Figure 4. - Concluded. Deflection of limiting 
streamline at surface with respect to potential- 
flow streamlines. 
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Mach number ahead of shock, 


(a) Cone half -angle, 20°. 


Figure 5. - Local skin-friction coefficient. 
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(b) Temperature ratio, T w /T, 1. 

Figure 5. - Concluded. Local skin-friction 
coefficient. 
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1 2 3 4 5 

Mach number ahead of shock, M^_ 

(a) Cone half -angle, 20°. 

Figure 6. - Lift on a yawed cone. 
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(b) Temperature ratio, T /t, 1. 

Figure 6. - Concluded. Lift on a yawed cone. 
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(a) Cone half -angle, 20°. 


Figure 7. - Pitching moment. 
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Mach number ahead of shock, M]_ 
(b) Temperature ratio, T W /T, 1. 


Figure 7. - Concluded. Pitching moment. 
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